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Abstract 

We consider an asymmetric version of a two-dimensional Coulomb gas, made 
up of two species of pointlike particles with positive +1 and negative —1/Q (Q = 
1,2, ... ) charges; Q = 1 corresponds to the symmetric two-component plasma and 
the limiting case Q — > oo is related to the one-component plasma. The system lives 
on the surface of a sphere, and it is studied in both canonical and grand-canonical 
ensembles. By combining the method of stereographic projection of the sphere onto 
an infinite plane with the technique of a renormalized Mayer series expansion it is 
explicitly shown that the finite-size expansions of the free energy and of the grand 
potential have the same universal term, independent of model's details. As a by- 
product, the collapse temperature and the Kosterlitz-Thouless transition point (in 
the limit of a vanishing hard-core attached to particles) are conjectured for any 
value of Q. 
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1 Introduction 



The universal finite-size properties of two-dimensional (2D) critical systems with short- 
range interactions among constituents are well understood within the principle of confor- 
mal invariance 0-0. For a finite 2D system of characteristic size R at critical point, the 
dimensionless grand potential (3VL [(3 = l/(kT) is the inverse temperature] has a large- R 
expansion of the form 

(3Q = AR 2 + BR-—\nR+ ... . (1.1) 
6 

The coefficients A and B of the bulk and surface parts are non-universal. The coefficient 
of the In i?-term is universal, dependent only on the conformal anomaly number c of the 
critical theory and on the Euler number \ of the manifold (with a smooth metric and a 
smooth boundary) to which the system is confined. In general, % = 2 — 2h — b, where h 
is the number of handles and b the number of boundaries of the manifold (x = 2 for a 
sphere, x — 1 f° r a disk, x = for an annulus or a torus). 



For lattice models and their continuum limits at critical point, the quantity ( |1.1|) is 
usually called the dimensionless free energy, with the notation /3F, but from a rigorous 
point of view those systems are usually treated in the grand-canonical ensemble (one sums 
over all spin configurations, particle occupations, etc.), with the grand potential as the 
generating functional. Here, we will consider a specific class of continuous fluid systems, 
which exhibit at any temperature the universal finite-size behaviour of type (|1 . 1| ) in both 
grand-canonical (fixed fugacity, the grand potential) and canonical (fixed particle number 
density, the Helmholtz free energy) ensembles. Although the two ensembles give the same 
thermodynamic behaviour in the thermodynamic limit, their finite-size corrections are 
expected to differ from one another. However, as will be documented in this work, at 
least for the sphere geometry of the confining manifold, the universal finite-size correction 
term is the same in both ensembles. 

The fluids under consideration are the 2D Coulomb plasmas of charged particles, in- 
teracting through a long-range logarithmic potential. The one-component version (OCP) 
of equally, say unit, charged point particles in a uniform neutralizing background is for- 
mally related to the fractional quantum Hall effect ||. The model is exactly solvable 



at the inverse temperature (5 — 2 by mapping onto free fermions. The symmetric 
two-component version (TCP) of oppositely ±1 charged pointlike particles undergoes the 
collapse of positive- negative pairs of charges at (3 = 2, which corresponds to the exactly 
solvable free-fermion point of the equivalent Thirring model [HI ~ I P- 0|| - If a small hard core 
is attached to the particles in order to prevent the collapse, the 2D TCP (or Coulomb gas) 



undergoes the famous Kosterlitz-Thouless (K-T) transition [|TTJ from a high-temperature 
conducting phase to a low-temperature insulating phase at around j3 = 4. 

The long-range tail of the Coulomb potential induces screening. Due to the screening 
effect, the bulk charge-charge correlations are short-ranged, exponential or even Gaussian. 
On the other hand, the same screening effect causes that the induced electrical-field cor- 
relations are long-ranged |T2|, jT3]. As a consequence, the free energy of the 2D OCP and 
the grand potential of the symmetric 2D TCP are supposed to exhibit, at any temperature 
of the conducting regime, a universal finite-size correction of type (|1.1|) ; this correction is 
opposite to the one occurring in the massless Gaussian field theory ||) which has c = 1. 
The arguments for a critical-like behaviour were first given for Coulomb gases with peri- 



odic boundary conditions []14}| , then for Coulomb gases confined to a domain by plain hard 
walls ||15||, by ideal-conductor boundaries [O], and finally by ideal-dielectric boundaries 



17] , [ 18 1 - In all cases, the explicit checks of the predicted finite-size universal behaviour 



were done at the exactly solvable inverse temperature (3 = 2 for both OCP and TCP with 
various geometries of confining domains. 

The general derivations of the universal finite-size correction term in refs. ||14||- [|18[ 



were based on plausible, but not rigorously justified, heuristic arguments. Very recently, 
a direct derivation was made for the specific case of Coulomb systems living on the surface 
of a sphere [where, since there is no boundary, B = in (|1.1|) 1. By combining the method 
of stereographic projection of the sphere onto an infinite plane with linear response theory 
(TCP, ref. []T9|) or with density functional approach (OCP, ref. |^0[|), the prefactor to the 
universal In R correction term was related to the second moment of the short-range part 
of the density-density direct correlation function csr for the planar system. This second 
moment, representing a new sum rule for the corresponding density correlation function, 



was evaluated explicitly for both OCP |21j and TCP |22| by using a renormalized Mayer 
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expansion [2^], p4[| , namely by observing a remarkable "cancelation property" of specific 



families of renormalized C5#-diagrams. The obtained results confirm quantitatively the 
prediction of conformal invariance ( |1 . 1|) for a critical system as we had c = — 1. 

In the present paper, we study a more general asymmetric 2D Coulomb gas made up of 
two species of pointlike particles with positive and negative charges of arbitrary strengths 
q a (a — 1, 2). Without any loss of generality, we can choose 

qi = +1 and q 2 = -1/Q (1.2) 

with Q = 1,2,... a positive integer. The model interpolates between the symmetric 
TCP (Q = 1) and the OCP, obtained as the limiting Q — > oo case after subtracting the 
kinetic energy of the 2-species. For the asymmetric Coulomb gas confined to a sphere, 
by combining the method of stereographic projection [|l^], |20j with the technique of the 
renormalized Mayer expansion |2~lf-||24j| we prove that the prefactor to the universal mid- 



term in ( |1 . 1| ) does not depend on the parameter Q. The same universal term is obtained in 
both canonical (free energy) and grand-canonical (grand potential) ensembles, although 
in the previous studies of Coulomb systems the choice of the ensemble seemed to be 
relevant: the exact solution at (3 = 2 and finite-size corrections were obtained strictly in 
the canonical ensemble for the OCP and in the grand-canonical ensemble for the TCP. 

From a technical point of view, the asymmetry in the strength of particle charges 
induces an asymmetry in particle correlations: two independent two-particle combinations 

11 = 22 and 12 = 21 in the symmetric TCP are replaced by the three ones 11, 22 and 

12 = 21. The method of stereographic projection |TjJ, |2(J is shown to imply a new 



"missing" interrelation among the second moments of species pair correlations. This, 
together with a generalization of the "cancelation" phenomena in the renormalized Mayer 
expansion, helps us to extend new sum rules for the plane OCP |21] and symmetric TCP 



22] to the general asymmetric plasma and thus to find the prefactor to the universal 



finite-size correction for the latter model. 



In their pioneering work about the asymmetric 2D Coulomb gas ||25|| , the authors 
predicted a successive series of temperatures at which partial recombinations (or col- 
lapses) of positive-negative pairs of charges occur, with a consequent impact on the short- 
distance behaviour of like-species correlations and on the equation of state (discontinuities 
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of dp/ dp, p is the pressure, were suggested). However, the exact solutions of thermody- 
namics for the symmetric Q = 1 |26|] and the first asymmetric Q = 2 [27] cases give only 
one singularity for the pressure at one collapse temperature, related to a short- distance 
instability of the configuration integral for a neutral cluster formed of one particle of 
species 1 and Q particles of species 2. The collapse temperature is explicitly derived for 
any Q in this work. If one attach a hard-core a to the particles, the system remains stable 
also below the collapse temperature, and can reach the K-T phase transition. In the limit 
of vanishing a, the presence of the K-T transition in the symmetric TCP is indicated by 



a divergence of the second moment of truncated density- density correlations ||19|| . Having 
at disposal the generalization of this second moment sum rule, we suggest the position of 
the K-T transition temperature for any value of Q. 

The paper is organized as follows. The method of stereographic projection in the 



grand-canonical |19j and canonical |20| ensembles is briefly generalized to the asymmetric 
2D Coulomb gas in sections 2 and 3, respectively. Section 4 deals with the renormalized 
Mayer expansion. The final results, conjectures and some concluding remarks are given 
in section 5. 



2 Grand-canonical ensemble 

On the surface of a sphere of radius R, two particles i of species <Tj and j of species aj 
interact via the Coulomb interaction — q ai q aj ln[(2i?/L) sin(^/2)], where L is a length 
constant set for simplicity to unity and 9ij is the angular distance of particles [ESJ. The 
particle charges belong to the two-component set given by equation (|1.2j ). Under the 
neutrality condition, the grand partition function reads 



I lg <j)=1 (2gsinj^ 
n^sin^ 



/3/Q 2 

x / dnf dnf . . . <> i%^f^ \> +..., (2.i) 



where z a (a = 1,2) are the species fugacities and dfi- cr ' ) is an angle element around the 
position of particle i of type a. The term in (|2.1|) containing N particles of type 1 and QN 
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particles of type 2 equals [ziz§ R 2 ^ X+ ®^ 1 *Q^] N times a dimensionless integral dependent 
only on f3 and Q. Thus, InS obeys the homogeneity relations 

^ Q 2 <9z 2 2(1 + Q) \ 4QJ OR K ' 

Let us denote by the number of particles of type a on the sphere of radius R; 

]\f{S) = ^jcR 2 x the constant density . Since 

= z CT ^-lnS, a = 1,2 (2.3) 



one gets the relations 



satisfying the electroneutrality condition J2a Qa n ^ = 0. 

For fixed fugacities z a and in the large- R limit, InS behaves like 

ln 3 = 0p (4nR 2 ) + f(R) + ... , (2.5) 

where p is the bulk pressure of an infinite (plane) system and f(R) is the leading finite-size 
correction. Substituting ( |2.5| ) into ( |2.4| ) one obtains, in the limit R — > oo, 

where n a is the number density of species a for the plane system. The equation of state 
follows immediately, 

»=»(i-£). ("> 

where n = 5Z<r n a is the total number density. The finite-size corrections to the species 
number densities then take the form 

- SoTQjf 1 -^)" 1 ^- < 2 ' 8b > 



Under stereographic projection of the sphere from its north pole onto an infinite plane 
tangent to its south pole 0, an area element R 2 dfl on the sphere is mapped onto d 2 r on 
the plane according to 

R 2 dQ = dV , (2.9) 

and the angular distance 9ij between two points on the sphere is expressible through the 
distance |t\ — rj\ of their projections on the plane by 

m sin t = 7 — Jy 1 r2 xi/ 2 • ( 2 - io ) 



J 



After the stereographic transformation, Q2.1| ) becomes the grand partition function of 
the underlying asymmetric Coulomb gas formulated on an infinite plane: besides the 
two-body interactions —q Gi (la 3 In |r» — rj\, there is a species-dependent one-body potential 
acting on each particle, 

m r) = 2 (i-£) ln (i + ^) j ^( 1 -£)ji, (,lla) 

«') - ^ (l - In (l + ^) ^ (l - ^ , (2.nb) 

The density shift of the planar system at the origin 0, due to the one-body potentials 
V a (r), is given by linear response theory, 

6n a (0) = - E/ & 2 r(h a {0)h a ,{r)f{3V al {r) , (2.12) 

where h a (r) = J2i ^a,ai^( r ~ r i) is the microscopic number density of species a at r and the 
truncated statistical average is taken in the unperturbed planar system. Let us introduce 
the notation 

I aa , = J d 2 r r 2 (n a (0)n a ,(r)) T (2.13) 

for the second moment of the species-density correlation functions, with the obvious sym- 
metry I aa i = I a i a . Since the density on the sphere equals to its projection on the plane 
at the south pole 0, it holds 

Sn a (0) = njf 5 - n a . (2.14) 



Inserting ( |2.11 ) into ( 2.12 ) and considering then ( |2.14|) in (|2.8|) , one finally gets 



f(R) = -ClnR, 

c = 2.(1 + 3) (i-^) 



! 12 + h h \ 



Simultaneously, as a by-product, an important interrelation among J's results: 



)-!){-¥■'. 



JL 

4Q 2 J \Q 



-^1 



22 — 1 12 



(2.15a) 



(2.15b) 



(2.16) 



3 Canonical ensemble 

In the canonical format, the stereographic projection maps the homogeneous particle 
densities on the sphere onto the inhomogeneous ones on the plane, 



n 



{S) 



n n \r) 



1 + _r?_V R- 
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(S) 



2R 2 



(3.1) 



The local chemical potential of particles a on the plane is given by 

where F is the free energy functional. According to the elementary fluid theory [^9 

5fi a (r) 5 2 f3F[n] 



(3.2) 



ft 



Sn ff >(r / ) 



5n a (r)5n a ,{r / ) 

-cw(r, r) + S a ^ , 



(3.3) 



where c stands for the direct correlation function. The chemical potentials of the asym- 
metric Coulomb system on the sphere coincide with the ones at the origin of the plane 
system, /v(0), for the density profiles (|3.1|). 

Let us start from a planar asymmetric Coulomb gas with particle densities and 
switch on the one-body potentials V a (r) (|2.11|) . According to ([3.1|) , the density change at 
point r caused by these one-body potentials is 5n a (r) ~ —n^r 2 /(2R 2 ). Using (|3.3[), the 
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induced shift in the chemical potential at the origin is 



(35^(0) = E/d 



2 r 



n a {0) 



8n a >(r) 



- ^E"?/^ c-Ary. (3.4) 

cr' 

Consequently, 

(3j2N^5fi a (0) ~ 27rJ2n { fH S > ) [ d 2 r c CT(7 ,(r)r 2 . (3.5) 

° Cr,CT' 

We will now omit the superscript (S) in the notation and denote N = ^Z a N a . The 
free energy F has the large- N behaviour 

(3F = AN + g{N) + const + . . . . (3.6) 

Since fi ff = dF/dN a , 

/mm. - ^ • (3.7) 



Taking into account that b\x a = 5fM a (0), eq. (|3.5|) then implies 



g(N) = DhiN , (3.8a) 
D = 27iJ2rian a/ J d 2 r c aa ,(r)r 2 , (3.8b) 

(j,<t' 

where n CT are the particle densities of the asymmetric Coulomb gas formulated on an 
infinite plane. 

4 Renormalized Mayer expansion 

We first review briefly the renormalized Mayer expansion in density (for details, see refs. 
|2T|j - PH ) for a general fluid composed of distinct species of particles {a} with the cor- 
responding densities {n(r,a)}. The particles % and j interact through the pair potential 
v(ri,<Ji\rj,<jj); vector position r*i will be sometimes represented simply by i. The renor- 
malization technique is based on an expansion of each Mayer function in the inverse 
temperature j3, and on a consequent series resummation of two-coordinated field circles; 
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by coordination of a circle we mean the number of bonds meeting at this circle. The 

renormalized -fT-bonds are given by 
K -(3v 



0~ww-0 = O- — -o + o-- 

1,01 2,0 2 1,01 2,0 2 1,01 



- -o + o- - 
2,02 l,0i 



- -o + . . . 

2,02 



or, algebraically, 



K(l,ai\2,a 2 ) = [-(3v(l, a x \2, a 2 )\ 

+ E /[-/Ml, 0i |3, 03)] n(3,0 3 ) ir(3,0 3 |2,0 2 ) d3 . (4.1) 

It is straightforward to verify by variation of ( }4.1| ) that it holds 

<y#(i,<7i|2,<7 2 ) 



X(l,0 1 |3,0 3 )X(3,0 3 |2,0 2 ) . 



(4.2) 



5n(3,0 3 ) 

The procedure of bond-renormalization transforms the ordinary Mayer representation of 
the excess Helmholtz free energy, defined by — /3F ex = — flF + Y^ a J dl [n(l, 0) lnn(l, 0) — 



n(l, 0)], into 



(3F e 



+ D {0) [n] +Y j D ( - s) {n] 



where 



/ \ 11 

> + / X + I I + 

* ^ - - .4 



is the sum of all unrenormalized ring diagrams and 



(4.3a) 



(4.3b) 



e 






etc. 



(4.3c) 



£)(1) £>(2) £)(3) £)(4) £)(5) 

represent the set of all remaining completely renormalized graphs. The free energy is the 
generating functional for the direct correlation function c in the sense that 



c(l,0i 2,0 2 ) 



Sn(l, 0i)<5n(2, 2 ) 



(4.4) 



With regard to QPD 



l,0i|2,0 2 ) = -^(l,0!|2,0 2 ) + -K 2 (l,0!|2,0 2 ) +^ C ^(l,0i|2,0 2 ) 



8=1 



(4.5a) 
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where c^(l, <ti\2, cr 2 ) represents the whole family of (1, cri)(2, <j 2 )-rooted diagrams gener- 
ated from a given D^ s \ 

p £)(«) 

C W(l,o-i|2,a 2 ) = — - — — r . (4.5b) 

dn(l, ai)dn(2, a 2 ) 

To get explicitly c^ s \ one has to take into account that the functional derivative of D^ s ' 
with respect to the density field generates root circles not only at obvious field-circle 
positions, but also on K bonds according to formula (|4.2|) , causing their right K-K 
division. The pair distribution function h, defined by 

n(l,a 1 )n{2,<r 2 )h{l,a 1 \2,a 2 ) = (n (T1 (l)n (72 (2)) T - 5^8(1 - 2)n(l, <n) , 

(4.6) 

is related to c via the Ornstein-Zernike (OZ) equation 

/l(l,<7i|2,<7 2 ) = C(l,<7i|2,<7 2 ) 

+ ]T / d3 c(l, <n|3, <r 3 )n(3, <r 3 )/i(3, <7 3 |2,<t 2 ) • (4.7) 

Let us return to the 2D asymmetric Coulomb gas with two kinds of particles a = 1,2 
of charges ( |1.2| ), interacting via the logarithmic interaction 

v(i,<Ti\j,<Tj) = q 0i q Vj v(i,j) , (4.8a) 
v(i,j) = —la\i — j\. (4.8b) 

We consider the infinite-volume limit, characterized by homogeneous densities n(l, a) = 
n a constrained by the neutrality condition 1cT n o = 0, so that 

rii = n — , n 2 = n — , (4.9) 

1 + Q ' l + Q K J 

with n being the total particle density. Two-body functions are both isotropic and trans- 
lationally invariant, c(l, <7i\2, <r 2 ) = c CTlCT2 (|l — 2|), h(l, o"i|2, cr 2 ) = h aia2 (\l — 2|). From eq. 
Q4.1|) it follows that the renormalized bonds exhibit the same charge-dependence as the 
interaction under consideration 



K{l,a 1 \2,a 2 ) = q ai q a2 K{l,2) , (4.10a) 
K(l,2) = -/3K (k\1-2\) . (4.10b) 
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Here, Kq is the modified Bessel function of second kind and k is the inverse Debye length 
defined by 

k 2 = 2vr /3 'E<?X • (4-11) 

CT 

The direct correlation function c, given by ( |4.5a|) , reads 

1 oo 

c aia2 (\l - 2|) = M(T2 /31n |1 - 2| + [^ (k|1 - 2|)] 2 + £ C ^ 2 (|l - 2|) 

Z - 8=1 



(4.12) 



For renormalized bonds of type ( fOO) , the c^-families of diagrams exhibit, as separate 
units, remarkable "cancelation properties", regardless of the topology of the generating 
graph D( s \ Let the given completely renormalized diagram (s = 2, 1, . . . ) be com- 
posed of U skeleton vertices % — 1, . . . , U of bond coordinations {z/j} and V = Y%=i v %l^ 
bonds a = 1, . . . , V. Omitting topological factors, can be formally expressed as 



£>(») 



u 

n 

i=i 



- v 

n^(«i 

. a=l 



(4.13) 



where a 1; a 2 G {1, . . . , £/}, a x < a 2 , denote the ordered pair of vertices joint by the a-bond. 
For the family & s \ generated from according to ( [4.5b ), one gets 



EwiJ(|r-^|) = n(E?> ff! 



u 



i=l \ <Ji 



(a) 

(b) 
(c) 
(d) 

(e) 

(/) 



x 



i i,j a 

+n | ndzE^-oE^Krwo^n^ 

i i a, b^a 

+" J II di E ^ " £ r)#(r, a 2 ) J] A b 

i i a- by^a 

j n dz e ^ k o ^ #(*■, «2) n ^ 
/ n di e ^ k r ) a 2 ) n k 6 
[ndi e ^(01,^^,02)^(61,0^,62) n ^ 



+ (n) 2 



(n? 



a,b 
(a#b) 



c^£a,b 



(4.14) 



where we have applied the functional relation ( |4.2|) . Here, K a = if(ai,a 2 ) defined by 
( |4. 10a| ) and n = J2aQa n ^ = n /Q- The (a) term in ([4.14j) corresponds to the creation of 
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root points at two skeleton vertices, the next (b, c) terms to one root circle generated at 
the skeleton and the other one at a bond, the (d, e) terms to two root points at the same 
bond and the last (/) term represents root points generated at two different renormalized 
bonds. Comparing (|4.14 ) with formula (31) of ref. [21] and assuming that k 2 = 27r/3n, it 



becomes clear that the rhs of (|4.14) is proportional to c^(\r — r 7 ]) generated from 



for the case of the 2D OCP with mobile-particle density h. As was proven in ref. |2T 
the zeroth and second moments of such a vanish, regardless of the topology of the 
generating diagram D^ s \ Consequently, in our model, 

5> CT rv/d 2 r cS(r) = 0, (4.15a) 

<7,cr' 

5> CT rv/d 2 r r 2 c%,(r) = 0, (4.15b) 

<J,(T' 

for every s = 1, 2, . . . . 

The 2D Fourier transform is defined by 

f(r) = -L fd 2 k /(fc)exp(ifc. r) , (4.16a) 

Z7T J 

f(k) = — [d 2 r /(r)exp(-ifc- r) . (4.16b) 

27T J 

For a function with circular symmetry, whose Fourier component is analytic in k around 
k = 0, one has 

oo 

/(*) = £/ (2i) ^\ ( 4 - 17a ) 

j=0 



/W = ^— / d 2 r r 2l f(r) . (4.17b) 
J 4^(j!) 2 27r ' 



The Fourier component of Kq(t) is given by 
d 2 r 



/(XT f°° 
Kq (r) exp (— ik ■ r) = / dr rJ (A;r)K 2 (r) 
27T Jo 

In (k/2 + y/l + (fc/2) : 



(4.18) 
fc^/l + (A;/2) 2 

with J being the ordinary Bessel function. The Fourier transform of the direct correlation 
function c, eq. (4.12), thus reads 



oo 



8=1 

(4.19) 
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where c^t(k) are analytic functions of k. Regarding the neutrality condition, relations 
( fOj| ) can be reflected via 



Y.n^Ak) = - ^ + O(k-) (4.20) 

Note that the coefficient to the /c 2 -term is universal, i.e., independent of (3, n\ and n-i- In 
fc-space, the OZ equation (4.7) takes the form 

Ka'(k) = c aa \k) + 2ir^c au n{k)n a „h a , la ,(k) (4.21) 



The insertion of the leading term c ac7 i(k) ~ — q a <ia'l3/k 2 into ( |4.21| ) leads to the zeroth- 
moment (electroneutrality) 



E^4t = -g ) *=1,2 (4.22) 



(0) Qa 

a' 

and the second-moment 

E qan a q a/ n a >h®, = p ( 4 - 23 ) 

Stillinger-Lovett sum rules fl32fl for the truncated charge-charge correlation function 
(p(0)p(r)) T , where p(r) = Y^a Q^a-( r ) is the total microscopic charge density at r [note 
the definition of h in ( ]Oj )]. 

To make use of ( |4.20| ), we introduce an auxilliary function 



U{k) = Y,n a ,c a , a {k) (4.24) 

cr' 

satisfying the relations 

= (4.25a) 

nj[ 2) +n 2 ff ] = -J-j. (4.25b) 

At the lowest order, the OZ equation ( |4.21| ) gives 



E = /?> + 2vr E ffnoXX ■ (4.26) 
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Combining O) with flOg) , /f 5 and / 2 (0) occur only in the combination (|4.25a|) , so they 
can be eliminated. Consequently, one gets explicitly 

l l + Q 



1(0) 



1(0) 
' i 22 



1(0) 
"■12 



2tto[1 - /3/(4Q)] 27m 
l l + Q 

27m[l - /3/(4Q)] 2vmQ 
1 



(4.27a) 
(4.27b) 
(4.27c) 



2vm[l - P/(4Q)] ' 

Without going into details, these formulae imply the compressibility sum rule 
i.e., the zeroth moment of the truncated density- density correlation function (n(0)n(r)) T , 
where n(r) = J2a n a{ r ) is the total microscopic number density at r. 
At the next k 2 order, the OZ equation ( (4.21|) gives 



E n^l = ti 2) + 2vr E fS^X + 2. E fS^h™ 



( 2 )„ 7,(0) 



(4.28) 



Inserting (|427D into |L28|), /f } and / 2 (2) occur only in the combination fl4.25b| ), so they 
can be eliminated. One obtains 



h { 3 + Qh% 



(2) 



h-12 + Q^22 



1 



l + Q 



9 ^-, im ^^ +Q1 ^- (4 - 29a) 



l + Q 



2vr(A (0) M 2 2 ) + Q/ 2 (0) /,g ) ) . (4.29b) 



96tt% 2 1 - P/(4Q) 

From eqs. (|4.6| ) and ( |4.17a| ) it follows that the second moments I aa >, defined by ( |2.13| ), 
are related to h^, in the following way 



(4.30) 



From the complete set of equations ( |4.29| ), ( |4.25a| ), ( ^.23| ) and the one ( |2.16| ) obtained by 
the stereographic projection, one finally finds that f± = (3/(8tt), fy ' = (3/{87iQ 2 ) and 



'-li 



'22 



'12 



2(3/?-8Q 2 )Q3-6Q 2 ) 
3/3tt(/3 - 4Q) 2 (1 + Q) 2 ' 

2Q 4 (3/?-8)(/3-6) 
3/3tt(/3 - 4Q) 2 (1 + Q) 2 ' 
2Q 2 [3p 2 - 2/3(6 - Q + 6Q 2 ) + 48Q 2 ] 
3/3tt(/3-4Q) 2 (1 + Q) 2 



(4.31a) 
(4.31b) 
(4.31c) 
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With regard to the definition of J's ( |2.13 ), the second moment of the number density 
correlation function is 

]2 2w„w^t i - 3/?(g - i) 2 /(2g 2 ) 



J d 2 r r 2 {h(0)h(r)) r - 



12tt[1 - /?/(4Q)P (4 ' 32) 
This is the generalization of the new sum rule derived for the symmetric Q = 1 TCP in 
refs. H and pi . 



5 Final results and conclusion 



Inserting ( [4.31| ) into ( |2.15b| ), one gets C = 1/3. With regard to ( p75|) and ( [2.15a| ), the 
dimensionless grand potential, f3VL = — In H, has the large- R expansion 

1 



f3Q = -I3p [AttR 2 ) + -\nR + const + . . . (5.1) 

This expansion has the universal lni? term like in ( |1 . 1|) . Taking the sphere value of x — 2, 
its prefactor corresponds to the conformal anomaly number c = — 1 as was expected. 
With regard to (|4.17a|) , it follows from ( |4.20|) that 



Y,n a n a > J d 2 r c aa \r)r 2 = — . (5.2) 



Consequently, D = 1/6 in (|3.8|), and (|3.6|) takes the form 



/?F = AA^ + -lnA^ + const + ... (5.3) 
6 

Since N = 4nR 2 n, the dimensionless free energy f3F has the same universal correction 
term as f3VL. Relations (|5.1| ) and (|5.3| ), valid for any asymmetric 2D Coulomb gas, are the 
main results of the present work. 

As was promised in the Introduction, we now derive an expression for the collapse 
temperature as a function of Q. On the base of the exact results for the lowest cases 
Q = 1 |^6| and Q = 2 ^7|, let us study the short-distance stability of the configuration 
integral of a neutral cluster formed of one particle of species 1, placed at the origin 0, and 
Q particles of species 2: 



nd V 1(i<j) = lln J ' . (5.4) 
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Here, A is a screening length of the Coulomb system and A is a short-distance cutoff. 
After the change of variables 7*j = Ar^, I\ is expressible as follows 

In the limit A — > 0, 1^ diverges as soon as 2Q — (3 /2 — (3 / (2Q) < 0, so that the asymmetric 
Coulomb gas is stable against collapse for [3 < /3 co \, with 

«- = ^tq ' (5 6) 

At the collapse point, although the density of particles (for a given fugacity) diverges, 
the truncated density correlation function (h( y O)h(r)) T is finite for r ^ 0. It was conjec- 
tured in ref. and subsequently in ref. |19| that the formula for the second moment of 
(n(0)n(r)) T , having the dimension of (length) , might be analytically extended beyond 
the collapse point. This second moment remains finite up to the K-T phase transition (in 
the limit of vanishing hard-core): its divergence at the K-T point is probably associated 
with a power-law decay of (n(0)n(r)) T at asymptotically large distances. With regard to 
the sum rule ( |4.32|) , we therefore suggest that 

/?kt = 4Q . (5.7) 

Note that [3 co \ < /3kt> as it should be. Increasing Q, both f3 co \ and /?kt tend to infinity, in 
full agreement with the fact that there exist neither the collapse phenomena nor the K-T 
phase transition for the one-component plasma. 

In conclusion, we have studied an asymmetric two-dimensional Coulomb gas confined 
to the sphere. In both canonical and grand-canonical ensembles, the finite-size expansions 
of the free energy and of the grand potential are explicitly shown to possess the same 
universal term, independent of the model's details. The equivalence of both ensembles 
smears out the previously indicated difference in the treatment of finite-size corrections for 
the one-component and two-component plasmas. Whether such an explicit calculation of 
universal finite-size corrections is possible also for other geometries of confining domains 
is an open problem. 
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